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Abstract 

We derive the effective energy density of thin membranes of liquid crystal elastomers as the 
T— limit of a widely used bulk model. These membranes can display fine-scale features both due 
to wrinkling that one expects in thin elastic membranes and due to oscillations in the nematic 
director that one expects in liquid crystal elastomers. We provide an explicit characterization 
of the effective energy density of membranes and the effective state of stress as a function of the 
planar deformation gradient. We also provide a characterization of the fine-scale features. We 
show the existence of four regimes: one where wrinkling and microstructure reduces the effective 
membrane energy and stress to zero, a second where wrinkling leads to uniaxial tension, a third 
where nematic oscillations lead to equi-biaxial tension and a fourth with no fine scale features 
and biaxial tension. Importantly, we find a region where one has shear strain but no shear stress 
and all the fine-scale features are in-plane with no wrinkling. 
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1 Introduction 

Liquid Crystal Elastomers are rubber-like solids that display unusual mechanical properties like 
soft elasticity and develop fine-scale microstructure under deformation. This material consists of 
cross-linked polymer chains where rigid rod-like elements (mesogens) are either incorporated into 
the main chain or are pendent from them. These mesogens have temperature-dependent interaction 
which results in phases of orientational and positional order Da 1331 ■ We refer to two phases: a 
high temperature isotropic phase, where thermal fluctuations thwart any attempt at order and a 
nematic phase, where the mesogens have a characteristic orientation but no positional order. This 
average orientation of the mesogens in the nematic phase is represented by a director. 

Nematic-elastic coupling is a key feature of these materials |32l [22]. The isotropic to nematic 
phase transformation is accompanied by a very significant distortion of the solid: typically elonga¬ 
tion along the director and contraction transverse to it. Further, the director can rotate relative to 
the polymer matrix. This novel mechanism induces a degeneracy in the low energy states associated 
with the entropic elasticity of the polymer network, whereby the material has a non-trivial set of 
nearly stress-free shape changing configurations. This degeneracy can lead to fine-scale microstruc¬ 
ture like stripe domains where the director alternates between two orientations in alternating stripes. 
Together, all of this gives rise to soft-elasticity m- 

A theory of nematic elastomers, and specifically the entropic elasticity associated with it, was 
formulated in Warner et al. [32], and was used to show the emergence of stripe domains and 
soft-elasticity. Mathematically, the energy functional is not weakly lower-semicontinuous resulting 
in possible non-existence of minimizers: briefly minimizing sequences develop rapid oscillations 
that result in lower energy than its weak limit. These rapid oscillations are interpreted as the 
fine-scale microstructure in the material. DeSimone and Dolzmann pHj computed the relaxation 
wherein the energy density is replaced with an effective energy density that accounts for all possible 
microstructures. The effective energy does indeed show soft elasticity, and can be used as by Conti 
et al. [9] to explain complex deformation patterns in clamped stretch experiments on nematic 
elastomer sheets Ell- 

Experiments on nematic elastomers, like the one highlighted, have largely been performed on 
thin sheets or membranes. These structures typically have instabilities such as wrinkling, and 
consequently membranes of usual elastic materials are unable to sustain compression and the state 
of stress is limited to uniaxial and biaxial tension. Thus elastic membranes have been described 
heuristically by theories like the tension field theory of Mansfield |25| , and such theories have been 
obtained systematically from three-dimensional theories [26. . 281 [23] . 

The goal of this work is to derive an effective theory of thin membranes of liquid crystal elas¬ 
tomers that accounts not only for the formation of fine-scale microstructure but also instabilities 
like wrinkling. An important insight that results from this is the possible states of stress in these 
materials. We find that like usual elastic membranes, membranes of liquid crystal elastomers are 
also incapable of sustaining compression, and the state of stress is limited to uniaxial and biaxial 
tension. Importantly, due to the ability of these materials to form microstructure, there is a large 
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range of deformation gradients involving unequal stretch where the state of stress is purely equi- 
biaxial. Consequently, a membrane of this material has zero shear stress even when subjected to a 
shear deformation within a certain range. 

We start with a three dimensional variational model of liquid crystal elastomers, derive the 
effective behavior of a membrane - a domain where one dimension is small compared to the other 
two - as the T—limit of a suitably normalized functional as the ratio of these dimension goes to 
zero following LeDret and Raoult [23j and others mm 110] . Our variational model is based on 
a Helmholtz free energy density that has two contributions. The first contribution captures the 
elasticity associated with the polymer matrix. Developed by Bladon, Terentjev and Warner [Ml E], 
it is a generalization of the classical neo-Hookean model to account for the local anisotropy due to 
the director. The second contribution following Frank [16] penalizes the spatial non-uniformity of 
directors, and has been widely used in the study of liquid crystals. In the context of liquid crystal 
elastomers, this penalizes domain walls - narrow regions that separate domains of uniform director. 
The competition between entropic elasticity and Frank elasticity, precisely the square-root of the 
ratio of the moduli k of the Frank elasticity to /i of the entropic elasticity - introduces a length- 
scale. It turns out (e.g. [33] ) that ~ 10 — 100 nm. Note that the thickness h of a realistic 

membrane is on the order of 1 — 100 ^m depending on the application. Thus, one has two small 
parameters, and one needs to study the joint limit as both \Jn/n and h go to zero, but at possibly 
different rates. We do so by setting k, = and studying the limit h —> 0. 

We find in Theorem |4.1| that the T—limit and thus the resulting theory is independent of the 
ratio Kh/h- This is similar of the result of Shu [29] in the context of membranes of materials 
undergoing martensitic phase transitions. In other words, the length-scale on which the material 
can form microstructure does not affect the membrane limit as long as it is small compared to the 
lateral extent of the membrane. Consequently, the T—limit we obtain coincides with the result of 
Conti and Dolzmann m who studied the case k = 0. In fact, our proof draws extensively from 
their work. Specifically, their result provides a lower bound and our recovery sequence is adapted 
from theirs. 

The T—limit is characterized by an energy per unit area that depends only on the tangential 
gradient of the deformation. It is obtained from the density of the entropic elasticity by minimizing 
out the normal component followed by relaxation or quasi-convexification. We compute this by 
obtaining upper and lower bounds, and provide an explicit formula in Theorem |5.1| (also shown 
schematically in Figure 1). It is characterized by four regions depending on the in-plane stretch: S 
a solid region where there is no relaxation, C a liquid region where wrinkling and microstructure 
formation drive the effective energy to zero, W a wrinkling region where wrinkling relaxes the 
energy and A4 a microstructure region where stripe domains relaxe the energy. The techniques 
employed here are in the same spirit as those employed by DeSimone and Dolzmann m in three 
dimensional nematic elastomers. 

We also study the oscillations related to the relaxation by characterizing the gradient Young 
measures associated with the minimizing sequences in Theorem 6.1 We show that the oscillations 
in the region A4 are necessarily planar oscillations of the nematic director and involve no out of 
plane deformation while those in the region W are characterized by uniform nematic director and 
wrinkling. 

We use the characterization of the gradient Young measure to define the effective state of stress, 


and show that this coincides with the derivative of the effective or relaxed energy in Theorem 7.1 


The Cauchy stress is given in (7.7): it is general biaxial tension in S, zero in £, uniaxial tension 
in W and equi-biaxial tension in A4. As described above, the unique attributes of liquid crystal 
elastomers give rise to this region of equi-biaxial tension compared to membranes of usual elastic 
materials. 
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This paper is organized in the following manner: In Section[2j we fix some notation and comment 
on background results which are used throughout the paper. In Section[3j we describe our model for 
nematic elastomers, a model which incorporates the entropic elasticity of the polymer matrix and 
an elastic penalty on the spatial gradient of the director. In Section [4j we derive our effective theory 
for nematic elastomer membranes based on a notion of T-convergence. In Section [5j we provide 
an explicit formula for the energy density in our effective theory. In Section [6] we characterize the 
microstructure in the aforementioned regions Ai and W. Finally, in Section [7j we conclude with a 
notion of stress in this effective theory and its physical implications. 


2 Preliminaries 


We gather here the notation and some background results which we use throughout the paper. We 
denote with W 1 the n dimensional Euclidian space endowed with the usual scalar product u-v = u T v 
and norm |ri| = V u T u. The unit sphere in is denoted by S n_1 and it is defined as the set of 
all vectors u E M n with |«| = 1. The space of m x n matrices with real entries is labeled with 
K rax “. When m > 1 we denote with 0(n ) the orthogonal group of the matrices F E M nxn for 
which FF 1 = F t F = /, where / is the identity in M nxn and with SO(n) the rotation group of the 
matrices F E O(n) with detE = 1. Letting now F E W nxn , adj s (E) stands for the matrix of all 
s x s minors of F, 2 < s < minjm, n}. In the case m = 3, n = 2, if 


/n 

fl2 \ I 

' (/ 21/32 - / 22 / 31 ) 

/21 

hi , then adj 2 (E) := 

— (/ 11/32 — / 12 / 31 ) 

/31 

hi J \ 

v (/ 11/22 — / 12 / 21 ) 


If u : M 2 —> M 3 is a smooth map, then adj 2 Vu is normal to the surface with equation {it(a;) : x E M 2 }. 
Letting F E M 3x2 , Q E 50(3), R E 0(2) by Proposition 5.66 QT] it follows that 

|adj 2 (E)| = |adj 2 (QE5)|. (2.1) 

Later in this paper we label the norm of adj 2 (E), with F a 3 X2 matrix, with 5 = 5(F) := \ adj 2 (E)|. 
Furtheremore, we simply write adj = adj 2 both when dealing with the adjugate of 3 x 3 and 3x2 
matrices. 

Finally, we state a version of the polar decomposition theorem: given any F E M 3x2 and any 
rectangular Cartesian basis, there exist Ai > A 2 > 0, Q E 50(3), R E 0(2) such that 

F = QDR (2.2) 

for 

Ai 
0 
0 

In fact, Ai, A 2 are the principle values of F. 



° \ 

A 2 ■ (2.3) 

0 J 


We now recall some concepts in calculus of variations (cf. |11]L We say that / : M mxn —> 
M U {—poo} is polyconvex if there exists a convex function g which depends on the vector M(F ) of 
all the minors of F such that f(F) = g(M(F )). In the case m = 3 , n = 2 then f(F) = g(F, adj(E)) 
with g : M 9 — y MU{+oo}. We say that / : M mxn —> Mu{+oo} is quasiconvex if, at every F E M mxn , 
we have 


/ f(F)dx < / f(F + Vu)dx 

0,l) n J( 0,l) n 


( 0 , 1 )" 
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for every u E W)’°°((0, l) n ,M m ). Note that the foregoing inequality holds for every D open and 
bounded subset of R n with \dD\ = 0 [3j- Finally, / : M mxn -> KU {+ 00 } is rank-one convex if 
t f(F + tR) is a convex function for all F, F E M mxn with rank(F) = 1. 

If a function / : M mxn — > RU {+ 00 } is not quasiconvex, we define f qc the quasiconvex envelope 
of / as 

f gc := sup{/i < f,h quasiconvex}. 

Analogously, we define f c , f pc , f rc as the convex, polyconvex and rank-one convex envelopes re¬ 
spectively of /. In the general case of extended-value functions, convexity implies polyconvexity 
and polyconvexity implies both rank-one convexity and quasiconvexity, but quasiconvexity alone 
does not imply rank-one convexity. Therefore, if / : W nxn —> M U {+ 00 }, we have 

jpe < jqc^ jpc < jrc (2.4) 

On the other hand, in the case of a real-valued functions, quasiconvexity implies rank-one convexity 
and hence, if / : M mxn —y M, we have 

f c < f pc < f qc < f rc ■ (2.5) 


We give an alternative representation formula for the rank-one convex envelope of a function / : 
l raxn —> M U {+00} 

K K 

f rc (F) := inf 1J2 W( F i) ■ E XiFi = Fi ) satisfy^} 

i i 


with A* > 0 and A* = 1. Family (A i,Fi) satisfies a compatibility condition here labelled with 
Hk and defined in DU Sec. 5.2.5]. In the same spirit we define semiconvex hulls of a compact set 
1C C M mxn . The set 


KP C = \F E 


r x " : f(F) < sup f{X) for all / : 

xeK. 


polyconvex | 


is the polyconvex hull of 1C. The quasiconvex hull JC qc and the rank-one convex hull K7 C are defined 
analogously. The lamination convex hull IC lc of /C is defined 

K lc =(fE M mxn : /(F) < sup f(X) for all / : R mxn ->1U {+ 00 } rank-one convex). 
Equivalently, lC lc can be defined by succesively adding rank-one segments (see M), i-e- 

OO 

K} c = U /c w 

i =0 

where JC° = JC and 


/C (i+1) = /C (i) U {F = AFi + (1 - A)F 2 : Fi,F 2 E F (i) , rank(Fi - F 2 ) < 1, A E [0,1]}. 
The relations between the different notions of convexity imply the inclusions (see [13] ) 

X lc C JC rc C K qc C JC pc . 


We refer the interested reader to mi and [23] for a discussion of all the different notions of convexity 
and their relations. 
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Finally, we introduce the notion of a gradient Young measure that characterize the statistics 
of the fine-scale oscillations in the gradients weakly converging sequences (cf. [24] ). We define 
a homogenous H 1 gradient Young measure to be a probability measure that satisfies Jensen’s 
inequality for every quasiconvex function / : M 3x2 —y M whose norm can be bounded by a quadratic 
function. Let M denote the space of signed Radon measures on M 3x2 with the finite mass paring 

</b/>= [ f(G)dn(G). 

J R 3 X2 

Then the space of homogenous H 1 gradient Young measures is given by 

M qc : = {v E M : |M| = 1, {u, f) > f((u, id)) 

V / : M 3x2 —y M quasiconvex with |/(G)| < C(\G\ 2 + 1)}. (2.6) 


3 Model of nematic elastomers 

Consider a nematic elastomer occupying a region 0 in its reference configuration, and assume that 
it is in its stress-free isotropic state in this configuration. Let y : —> M 3 describe the deformation 

and n : 0 —> § 2 describe the director field. We denote Vrt to be the reference gradient of some field 
u : tl —> M 3 and V y u to be the spatial gradient of u. It follows V y u = Vu(F) -1 where F = Vy is 
the deformation gradient. 

We take the Helmholtz free energy density of the nematic elastomer to be the sum of two 
contributions: 

w = w e + w n 

where the W e describes the entropic elasticity of the underlying polymer chains of the nematic 
elastomer and W n describes the elasticity of the nematic mesogens. 

Following Bladon, Terentjev and Warner in |6, |33j, we take the entropic elasticity to be of the 
form 


f-fTr (F T r l 
W e (F,n) = { 2 V v 

I +oo 


F) -3 


if det F = 1, \n\ = 1, 
else, 


(3.1) 


where 


i = r 1 / 3 (I + (r — l)n <g) n) 


(3.2) 


is the step-length tensor. Here y > 0 is the shear modulus of rubber and r > 1 is the (non- 
dimensional) backbone anisotropy parameter. Note that for r = 1 the energy reduces to that of a 
neo-Hookean material. As most nematic rubbers are nearly incompressible |33j, we prescribe W e 
to be a finite only for volume-preserving deformations. We can substitute for i and write 


t f r l/3 


W e (F,n) = { 2 

+oo 


(l^l 2 


r — 1, 


IF 1 nr ) -3 


if det F = 1, |n| = 1, 
else. 


(3.3) 


For future use, we define a purely elastic energy by taking the infimum over directors. Following 
DeSimone and Dolzmann 


= '( de ‘ F = 1 ' 
«6 S 2 I +oo else, 


(3.4) 
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+oo 



where 




r — 1 


\2 

a m 


-3 . 


(3.5) 


Here A m{F) is the largest eigenvalue of (F T F) 1 / 2 . Energy density Wsd{F) is not quasiconvex and 
the quasiconvex envelope W^' D has been computed in [ 13 ]. 

Following Oseen, Zocker and Frank, (see for example, [15]). we take the elasticity of the nematic 
mesogens to be of the form 


W n = -fiq(divn) 2 + - K 2 (n • curln) 2 + -K%{n x curln) 2 


(3.6) 


where divn and curln are the spatial divergence and curl of the director respectively, and k\ > 
0, K 2 > 0, K 3 > 0 are known as the splay, twist and bend moduli respectively. Notice that this is a 
non-negative quadratic form in n and X7 y n. It turns out that these moduli are very close to each 
other and one can introduce an equal modulus approximation 


W?q = -\V y n\ 2 = -\(Vn)F 


—112 


(3.7) 


eq 2 ' ' y “' 2' 

where the second equality holds formally. Importantly, from a mathematical point of view, any 
given W n of the form (3.6) can be bounded from above and below by equal moduli approximations, 
and therefore all the results we prove for the equal modulus approximation hold for the more general 
form. Finally, since we assume incompressibility or det F = 1, F _1 = adj F so that 


W^ = -|(Vn)(adjF)| 2 . 


(3.8) 


Putting these together, the Helmholtz free energy of our nematic elastomer, F : i7 1 (H,M 3 ) x 
iy 1,1 (H, § 2 ) ->1U {—poo}, is given by 


F(y,n) := 


(3.9) 


J (V e (Vy,n) + ^|(Vn)(adj Vy)| 2 ) dx if (y,n) € A, 

_+oo else , 

where for definiteness, the set of admissible fields is 

A := {(y, n) G H\n,R 3 ) x W^faS 2 ) : (Vn)(adj Vy) G L 2 (H,M 3x3 )}. 

The ambient space for deformations, y G I7 1 (H,M 3 ), is optimal since for F G M 3x3 satisfying 
det F = 1, W e satisfies the growth and coercivity 


-\F\ 2 — c < W e (F, n) <c(|F| 2 + 1) 
c 


(3.10) 


independent of n G S 2 . Here c > 1 depends on r and fi. The ambient space for the director held, 
n G H /1,1 (H,S 2 ), may not be optimal. Nevertheless, consider the following: 

Remark 3.1 Fonseca and Gangbo m showed the lower-semicontinuity of 

[ \(\/n)(Vy)- 1 \ 2 dx 
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in the space 


Pp,q = {(u, y) £ W 1,p (fl, 8 2 ) x VF 1,9 (f2,K 3 ) : detVy(x) = 1 a.e. x E H}, 

where 2 < p < +oo and 4 < q < +oo such that (1/p) + (2 /q) = (1/2), (see Theorem, 4-1, W- The 
existence of minimizers in /3 Pi9 follows from this result. However, notice that this is more regularity 
than we assume. In fact, the existence of minimizers in A is not clear. However, this does not 
affect T—convergence or the membrane limit. 

Barchiesi and DeSimone showed well-posedness for an energy similar to \S. ,9[ j where n was 
taken as a mapping from the deformed configuration y(Sl), and it was assumed W e (F , n) > c(|F| 3 — 
1)- 


4 Membrane Theory 

In this section, we derive a theory for nematic elastomer membranes whose three dimensional free 


energy satisfies (3.9). 


4.1 Framework 

We consider a nematic elastomer membrane of small thickness h which has a flat stress-free isotropic 
reference configuration := { (x 1 , x 3 ) E M 3 : x' E u, x 3 E (—h/2, h/ 2)}. We assume co is a bounded 
Lipschitz domain in M 2 . Let y : fl/, —> M 3 describe the deformation and n : Dh —t S 2 describe the 


director field so that J~ h (y,h ) is the Helmholtz free energy in (3.9) now parameterized by the 
thickness of the membrane in its reference configuration. We assume n/2 = Kh, Kh h 0 and Kh —> 0 
as h —> 0. 

To take the limit as h —> 0, we follow the theory of T-convergence in a topological space endowed 
with the weak topology. The general theory can be found in [7] and [T2j. In order to deal with 
sequences on a fixed domain, we change variables via 

z = (zi, z 2 ) = (xi,x 2 ) = x', z 3 = -j-x 3 , x E Q h 

and set H := oj x (—1/2,1/2). To each deformation y : —> M 3 and director field h : fl/j —> S 2 , 

we associate respectively a deformation y : fl —> M 3 and director field n : H —> § 2 such that 


y(z(x )) = y(x) and n(z(x)) = h(x), x E flh- 
We set T h (y,n) := T h (y,n)/h , and following the change of variables above observe 

- rhfn ^ _ ) (w e (V h y,n) + ^|(Vn)(adj Vy)| 2 ^) dz if (y,n) E A, 


(4.1) 


I ( y,n ) = 


(4.2) 


+00 


else 


where V/j y = (V'y\(l/h)d 3 y) with V' the in-plane gradient. We also use the identity (Vft.n)(adj V^y) 
(l//i)(Vn)(adj Vy). 

Finally, we take our membrane theory to be the T-limit as h —> 0 of the functional defined on 

H\n,] 


i h (y) := 


inf 

neiWTr^s 2 ) 


l h (y,n). 


(4.3) 




4.2 The Membrane Limit 


Theorem 4.1 Let T h be as in (4-3) with Kh > 0 and Kh —> 0 as h -G 0. Then in the weak topology 


of T h is equicoercive and T-converges to 


J{y) = 


W 2 D^y'y) dz ' ifd 3 y = 0a.e., 


Too 


otherwise . 


Here 


W 2D {F) := inf W 3D (F\c) 

ceiR 3 


for W;io given in (3-4) and is the quasiconvex envelope ofW^D, 


(4.4) 


(4.5) 


Wl c D {F) = inf { [ W 2D {F + Vcftdz' : cf € ^’“((O, l) 2 , M 3 ) 1 . 


Equivalently: 

(i) for every sequence {yh} C 4t 1 (0,M 3 ) such that T h (yh ) < C < Too, there exists a y E 
H 1 { fl,M 3 ) independent of z 3 such that up to a subsequence 

yh ~ T Vhdz — 1 y in H l (n,R 3 y, 

Jo. 

(ii) for every {yjf] C such that yh —>■ y in i4 1 (J7,M 3 ), 

liminfX h (y /l ) > J(y); 
ft—H) 


(in) for any y G i4 1 (fl, 
iT 1 (ST, M 3 ) and 


t/iere exists a sequence {?//,} C i4 1 (12,M 3 ) sncd that y^ —*• y in 


limsup Z h \yh) <J{y). 
ft-> o 


The result for the case Kh = 0 was provded by Conti and Dolzmann [10] (Theorem 3.1 there). 
Theorem 4.2 (Conti and Dolzmann [10] ) In the weak topology o/14 1 (n,M 3 ) ; the functional 

Fe(y) ■■= [ W 3D (V hy )dz 


is equicoercive and T-converges to J given in ( 4-4 )■ 

Remark 4.3 A different dimension reduction theory for hyperelastic incompressible materials was 
developed by Trabelsi f3^,j3jy under similar assumptions. Trabelsi shows that the membrane energy 
density (integrand of J) is given by ((W 2 D) rc ) qc ■ From the proof of Theorem 5.1 below, it follows 
that Wff = W(f D (and hence {W.ff ) ) qc = W(f D ). Thus the two limits agree. 
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Remark 4.4 We remark on some general properties of the purely elastic portion of our nematic 
elastomer energy density. Wo : M 3x3 —> M in (3.5) is Lipschitz continuous, Wo + 3/i/2 is non¬ 
negative, and there exists a constant c such that 



The energy W 2 n in (4-5) is given by 


c < W 0 (F) < c{\F\ 2 + 1). 


W 2D (F) 


min W;$o(F\c) 

ceM 3 

+oo 


if rankF = 2, 
else , 


(4.6) 


(4.7) 


and satisfies 

\ (^ |2+ jtw) - c - W2D(P) - c i lP ' 2+ ih +1 ) m 

with 5(F) = |adj(F)|. The effective energy density W$p is quasiconvex, Lipschitz continuous on 
bounded sets and its definition does not depend on the choice of the domain u, as long as it is open, 
bounded and \du\ = 0. Furthermore, there exists (Lemma 3.1, f l(Jf ) a constant d such that 


^\F\ 2 -d <Wf D (F)<c'\F\ 2 + d. 


(4.9) 


Proof of Theorem f.l 


Note that trivially, 


l h (y) > [ inf W e (V h y, n)dz = lf(y). 
Jn ne§ 2 


(4.10) 


Therefore, the compactness and lower bound (Properties (i) and (ii) in Theorem 4.1) follow from 


Theorem 4.2 It remains to show Property (iii). This is done in Proposition 4.6 


Before we proceed, we note that the fact that the T—limit is independent of Hh/h is similar to 
the following result of Shu (29j . He also provides some heuristic insight. Since the membrane limit 
optimizes the energy density over the third column of the deformation gradient, there is little to 
be gained by oscillations parallel to the thickness. Consequently, penalizing these oscillations with 
Kh does not affect the T—limit. 


Theorem 4.5 (Shu |29| ) Let Kh —> 0 as h —> 0, and W : Ll —> M be continuous and bounded 
from above and below by \F\ P ± c respectively for some c. Then, in the weak topology of W 1,p , the 
functional 


{«h|V ft V A y| 2 + W(X7 h y)}dz, 


T—converges to 


L{ W 2o( V 'y)} dz ' l f d 3V = 0 a - e -> 

else , 


JUJ 

oo 


where W 2 d(F) := inf ceR 3 W 3D (F |c). 


10 










4.3 Construction of Recovery Sequence 

It remains to construct a recovery sequence to prove J is the T-limit to T h . 

Proposition 4.6 For every y £ i4 1 (f2,M 3 ) independent of Z 3 , there exists a sequence {( yh,nh )} C 
C°°(li,M 3 ) x C 1 (Q,S 2 ) such that yh —*■ y in i4 1 (H,M 3 ) and 

limsup T h (y h ,n h ) < J(y). (4.11) 

h-> 0 

Our construction also draws heavily from Conti and Dolzmann (lOj . The main difference is that 
we need additional regularity for our recovery sequence We summarize the Conti-Dolzmann 
construction in two lemmas. The first lemma regards the construction of a sequence to go from 
the energy density W 2 d to W,f D on For our analysis, the important observation is that in the 
limit the deformation gradient is constant on an increasingly large subset of oj. The second lemma 
regards the extension of smooth maps on u to incompressible deformations on f 

Lemma 4.7 (S. Conti and G. Dolzmann [10]) For any y £ H 1 ( w,M 3 ) ; there exists a sequence 
{yj} C C'°°(w,M 3 ) such that rank Vyj = 2 everywhere, yj —*■ y in Ft 1 (oj, M 3 ) as j —> 00 , and 

limsup f W 2 D (Vy j )dx' < f W^(Vy)dx'. (4.12) 

j — y 00 J to J to 

Moreover, the sequence has the following properties: 

(i) for each j £ N, yj is defined on a triangulation T - 7 of u which is the set of at most countably 
many disjoint open triangle T- J whose union up to a null set is equal to oj, and P,- is the jump 
set given by 

Tj := du U (J dTM, 

i 

(ii) there is a sequence of boundary layers {r]j} such that rjj > 0 and rjj —> 0 as j ^ 00 , and the 
set T rij is defined to be 

T Vj := {x' £ oj : dist(x / , Tj) < rjj}; 

(in) if TC \ T rij is nonempty, then V'yj is a constant on this set and we set 

F/ :=Vy,(*'), x'£T/\T % ; (4.13) 

(iv) adj V'yj is bounded away from zero in the sense that for some €j > 0 sufficiently small, the 
inequality 

I adj V'y-j | >ej> 0, 

holds everywhere. 

Lemma 4.8 (Conti and Dolzmann [TO]) Let w, v £ C 00 ^, M 3 ) satisfy 

det(V / rc|^) = 1 in w. 

Then there exists an ho > 0 and an extension v £ C°°(a) x (—ho, ho), M 3 ) such that v(x', 0) = w(x') 
and detVu = 1 everywhere. Moreover, for all X 3 £ (—ho, ho) the pointwise bound 

|Vv(x) — (V , tc|i/)(x / )| < C\xf\ 

holds, where C can depend on w and v. 
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We construct a recovery sequence and thereby prove Proposition 4.6 in 4 parts. In Part 1, we 


take a sequence of smooth maps yj as in Lemma 4.7 and show that we can construct a sequence 
of smooth vector fields Cj such that det(V / y J |cj) = 1 in uj. In Part 2, we use Lemma 4.8 to extend 
yj appropriately to a deformation on II, i.e yj \ In Part 3, we construct a sequence of C 1 director 
fields nj 1 on 12 which enables passage from W e to II '2 n ■ Finally, in Part 4 we show that we can 


take an appropriate diagonal sequence hj —> 0 as j —> 00 which proves Proposition 4.6 


Proof of Proposition f.6. Let y E H l (Q, 
(with abuse of notation). By Lemma 4.7 
rank V'yj = 2 everywhere, yj —*■ y in H\uj, 
and the sequence satisfies properties (i 


independent of z 3 . Then y is bounded in 77 1 (w, 


we find a sequence {yj} C C°°(uj, M 3 ) such that 


the energy is bounded in the sense of (4.12), 
(iv) from the lemma. 

Part 1. We define the smooth vector field Cj on the triangulation T - 7 for yj in Lemma 
each nonempty T 3 \ r, 7j . there exists a constant F 3 defined in Lemma 

Motivated by t 


4.7 


Then by (4.7), Wsd^FPc) has a minimizer for c E 


4.7 


(i). On 
(iii), and it is full rank. 


lis observation, we let 


which via (3.4) implies 


Consider the vector field, 


= arg min W 3 d ( F \ 3 c ), 

ceR 3 

(4.14) 

det(F/|c/) = 1 . 

(4.15) 

j adj Vyj 

!i ladjV'^l 2 ' 

(4.16) 


This is well-defined given Lemma 4.7 (iv). Moreover, since yj is smooth, Cq E C°°{ cj, 
since det(F|c) = (adj F) T c, we have 


Let Cj E C' 0 O (u;,M 3 ) be given by 


c o + fa c i ~ c i 


det(V‘ 2 /y |cq ) = 1 in uj. 

on each T/ \ T Vj with nonempty open subsets, 
otherwise on uj. 


Further, 


(4.17) 


(4.18) 


Here ipi E 6 ( 5 °(T )/ 7 \ T r]jl [0,1]) is a cutoff function which equals 1 at least on the entirety of the 
subset Tj - 7 \ r 2)?J . Notice when Cj = c 0 J , the determinant constraint is satisfied trivially by (4.17). 
Conversely, combining (4.15) and (4.17), 


det(Vyj\cj) = (adj Vyfy (c+ if, 


C- - Cn 


= det(V / y i |c ( j') + ifi (^det(F. J |c/) - det(V , y i |c 0 : ')^ 

= 1 on each T ri 3 \ T v with nonempty open subsets , 


since V'yj = F { J on this set. We then conclude detfV' yj\cj) = 1 in uj, and this completes Part 1. 
Part 2. Fix j E N. From Part 1 we have yj, Cj E C' 00 (d), M 3 ) satisfying det(V'y^jcj) = 1 in uj . He nce, 
there exists an > 0 and avE C°°(uj x h(})) such that the properties of Lemma 


4.8 


hold 
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replacing w with yj and u with Cj. Let h E (0, h 0 J ) and y- 1 E C°°(£lh, 


Qh- Further, let yj 1 E C 00 ^,!! 3 ) be associated to y,j' using (4.1). From Lemma 
y/(z',0) = yj(z'), det V h y^(z) = 1 and 

I VhyHz) ~ (V' yj\cj)(z')\ < Cjh\z 3 \ < Cjh, z&VL. 


be the restriction of v to 
we conclude 


4.8 


(4.19) 


Here Cj is a constant depending on yj and Cj. and the second inequality above follows since 
z 3 E (—1/2,1/2). From these properties we conclude as h —> 0, 


1 


yJ 1 —> yj in 1L 1 (0,M 3 ) and -d 3 yj 1 —> Cj in L 2 (H,M 3 ). 


h 

This concludes Part 2. 

Part 3. As in Part 2, we keep j E N fixed. From Lemma 


(4.20) 


4.7 


(i) we have that T 3 = w (up to a 
set of zero measure), though this union can be countably infinite. From herein, we choose a finite 
collection of N(j) triangles so that 


L S% < {wM + 1}dz ^T 

Then for each of the N(j) triangles for which the set T 3 \ T ri - is nonempty, let 

n- := arg min W e (F- 3 1 c- , n ). 

* nes 2 * 1 

Further, let q 3 be the piecewise constant function on M 2 given by 


Qj(z') ■= 


n/ if i E {1,..., N(j)}, T 3 \ F Vj is nonempty, and z 1 E , 
q otherwise in M 2 . 


(4.21) 


(4.22) 


(4.23) 


Here q is a fixed vector in § 2 . Then q 3 maps to S 2 , but it is not in C 1 . To correct this, we employ 
the approach used by DeSimone in m (see Assertion 1). 

Observe by construction the range of q 3 is finite. Hence, there exists an Sj E § 2 and a closed 
ball B e (sj) of radius e > 0 centered at Sj such that (range qj) 0 B e (sj) = 0. Then the stereographic 
projection Tr Sj with the projection point as Sj maps the range of qj to a bounded subset of K 2 . Let 

(ii), and consider the composition 


be a standard modifier with rjj as in Lemma 


4.7 


hj := tt 1 o 


Vi * t 71-5 i 0 


Qj))- 


This composition is well-defined since the range of qj is outside a neighborhood of the projection 
point Sj. Further, hj maps to 8 2 using the definition of the inverse of the stereographic projection. 
Moreover, 7 r” 1 is differentiable and its argument 7 / 7 ^- * (7r Sj o q 3 ) is smooth. Hence, hj E C' 1 (M 2 ,S 2 ). 

Let rij E C7 1 (a>, S 2 ) be the restriction of hj to the closure of ui. Further, let nj 1 E C' 1 (D,§ 2 ) be 
the extension of rij to H via nj l '(z) := nj(z') for each zEU. As a final remark for this part, observe 
for * E {1,..., N(j)} and z' E T 3 \ T 2rij , 

n h (z) = nj(z') = 7 r 


-1 

’3 


O * (n. 


= 7 T„ 1 O 


Vj * V /Is 7 u Qj)) ( z ) 

^Vj( z> - ° Qj)(£)d£ 


= n Si L O [ (7 T Sj o qj) / ^ (z 1 -£,) d Z 

J B V j( z ') j 

= TT^ 1 0 0 Qj) = n i, 


13 


(4.24) 






since B, r]j (z') n dTd = 0 and so qj is constant on B r/j (z > ), see (4.23). This completes Part 3. 

Part 4. From Parts 1-3, we have for each j E N the functions yj 1 E C 00 ^,!! 3 ) and nj 1 E C 1 (fl,S 2 ) 
parameterized by h E (0,/ig). It remains to bound the functional X h appropriately and take the 
lirnsup. For the bounding arguments, C shall refer to positive constant independent of h and j 
which may change from line to line. From (3.1), when W e is finite, it satisfies a Lipschitz condition 

| W e (F, n) — W e (G,n)\ < \£~ 1/2 | 2 (|F| + |G|) \F - G\ 

<C(\F\ + \G\)\F-G\. 

As asserted above, \l~ 1 d 2 \ is uniformly bounded for n E § 2 . Then since for every z E 12, det (V f l yj l )(z) 
1, det(V , y J |c :; )(z / ) = 1 and n/(z) = nj(z! ) E § 2 , 

[ W e (V h yj l ,nj l )dz < [ W^Vy^^dz' + [ \W*{V h yf, n/) - W e (V yj \ Cj , nf)\dz 
J fl Ju J 

< [ W e (V , y j \cj,rij)dz' + E^j. (4.25) 


Here E^j is the estimate obtained from the Lipschitz condition and an application of Holder’s 
inequality, 

E h,j '■= C (||V/»i^ A ||x,a ( n jR 3 ) + ||(V , yi|c :/ -)||L 2 (Q )K 3 )^ \\VhyJ 1 - (y'yj\cj)\\L 2 (ci,R 3 )- (4.26) 


We now focus on the first term in the upper bound (4.25). For i E {1,..., N(j)} and z' E 
Td \ r 2 V:I , observe 


W e (y' V j(z')\cj(z') , rij{z')) = \V (Fd'.cj 


by 


by (4.22 


= min W e (F /1 c/, n) 
ne § 2 4 1 1 ' 

= min W 3D (Fd\ c ) 
ceR 3 1 

= W 2D (Fd ) = W 2D {Vy j {z 1 )) by <|47h. 


4.13), (4.18) and (4.24); 


by (3.4) and (4.14); 


Then, 


W e (Vy j \c j ,n j )dz' < 


+ 


' . W 2D (V' yj )dz' + 

(ufjfTdm^ 

W e (V'yj\c j ,n j )dz', 


L\u N< -p T d 

' 4 = 1 X 


W e iy'yj | Cj , nj)dz' 


(4.27) 




using our result for W e and since each integrand is nonnegative. 

W e bound W e (V'yj\cj, rij) in (4.27). To obtain thi s bo und notice \c(j\ 2 = 1/| adj V / y ? -1 2 from 
(4.16). Further, using the coercivity condition of Wq in (4.6), the definition of cj in (4.14), and the 
growth in (4.8), 


c /| 2 < Wo{Fd\cd) = W 2D (Fd) < c I \Fd\ 2 + 


j 12 


adj F- 


J\2 


+ 1 
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Following these observations, we notice on the sets T- \ T^., X7yj = F - J by definition (see Lemma 


4.7 (iii)) and therefore, 


c j\ 2 < 2(|Cq | 2 + |c/| 2 ) < C ( \V'yj\ 2 + 


adj V'yd 2 


+ 1 


since Cj is as in (4.18). On the exceptional sets, by definition cj = c 0 , and the right side above is 
still an upper bound to \cj\ 2 . Hence everywhere in u, 

W e (| Cj , rij ) < c (|VVi| 2 + lof + 1) 

5 ° ( |V '» |! + adjV'fcf + ') 

<c(w 2D (V yj ) + 1 ), 


using the growth in (|3.10|), the bound above and the coercivity in (|4.8|). This implies the bound 

(4.28) 


W e {Vy j \cj,n j )dz' < / W 2 D(Vy j )dz' + E 2 , 


where recalling (4.27) and (4.21), the remainder E 2 is given by 


E],= C 


1 


/ {W 2 D {yy j ) + l}dz l + -) . 
' r 2n, 3 J 


(4.29) 


To recap, from (4.25) and (4.28), the entropic part of the energy is bounded above by the 
estimate 


W e (V h yf l ,n : h )< [ W 2D (Vy j )dz' + El j + E 2 . 
! J u 


(4.30) 


It remains to bound the elasticity of th e nem atic mesogens. 

Consider the second term of X h in (4.2). Our deformations and director fields have sufficient 
regularity, so 


I |(Vn/)(adj V?//)| 2 cfe = n h I |(V / n J |0)(adjV ft y/)| 2 cte. 


(4.31) 


Here, we used the identity (l/fi)(Vn)(adj Vy) = (V/ l n)(V/ ! ,y) and the definition n^z) := nj(z'). 
We bound the integrand by a constant independent of h. To do this, we first consider the pointwise 
estimate in (4.19). An application of the reverse triangle inequality on this bound yields for small 
h the pointwise estimate 


|Si»/(2)l 2 + l&y/WI 2 + p l*»/WI 2 = NkS^z)]- 


< (Cjh + l(V'*/,|c,)(/)|) 2 

< [Mj! 3) 1/2 , zen. 


(4.32) 


Here Mj is a constant which depends only on yj and Cj. Then F = (/ 1 I/ 2 I/ 3 ) £ M 3x3 satisfies 

I adj F\ 2 = |cofF| 2 = |(/ 2 x / 3 |/ 3 X /i|/i X h)\ 2 

= I /2 x / 3 | 2 + |/ 3 x fi\ 2 + \fi x f 2 \ 2 
<l/ 2 | 2 l/ 3 | 2 + |/3| 2 |/il 2 + l/il 2 l/2| 2 , 
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and so we can bound from above (4.31) 


Kh [ |(V , nj|0)(adj V h y)\ 2 dz < K h [ \Vnj\ 2 \ adj Vny^dz 
J Q J n 


< Kh l t iv'n^l 2 ^J^\d 2 yj l \ 2 \d 3 yj l \ 2 + ~^ 2 1 831 /j 1 1 2 1 diyj 1 | 2 + |3iZ/j| 2 |d2l/j| 2 ) dz. 


Applying the bound in (4.32) to this estimate, we conclude as desired 




/ |(adjVn/)(Vj//)| 2 d 2 < Kh Mj / \V rijfdz =: n h Mj 


(4.33) 


Here Mj is a constant depending only on y-j, Cj and rij. 

To complete the proof of Proposition 4.6 it remains to show that in the limit as h — > 0, the 
energy is bounded as in (4.11). From (4.30) and (4.33), 


< / W 2D {Vyj)dz' + El , + E 2 + n h Mj. 


(4.34) 


We now fix j e N and take the limit as h —»• 0. Notice from (4.20), \\S7hVj ~ (V / r/.,-|cj )|| i 2 0 

as h —> 0. This implies W^hV^WL 2 5- Cj for some constant Cj independent of h. With these 
two observations, we conclude £)C 


0 as h —^ 0, see (4.26). Further, since Kh —>• 0 as h —» 0, 


KhMj —> 0 since Mj is independent of h. Collecting these results and combining with (4.34), 
limsup i h \yK n- 1 ) < limsup ( f W 2 D(y'yj)dz' + E\ ■ + E 2 + KhMj) 

h —^0 h-+0 \Ju i J 

= [ W 2D (V yj )dz' + E 2 . 


Finally, using (4.12), the fact that |r 27 j-| —> 0 as j —> 00 (rjj -» 0, see Lemma 4.7 (ii), and (4.29)) 
we conclude 


limsuplimsup T h (yj 1 , nj 1 ) < limsup / W 2 n(S'yj)dz! + E 


j — ^00 h—t 0 


< [ W? D {Vy)dz'. 


We now choose a diagonal sequence hj —> 0 as j —> 00 so that this estimate is satisfied and y hj —*■ y 
in i7 1 (H,M 3 ). This completes the proof. ■ 


5 Effective energy of nematic elastomer membranes 

In this section, we provide an explicit formula for the effective energy of nematic elastomer mem¬ 
branes. 


5.1 Effective energy 


Theorem 5.1 Let W e as in (3.3). For any F e 


3x2 


, let 


W 2D (F) = inf W e ((F\c),n). 

cSK 3 ,neS 2 


(5.1) 
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Figure 1: The effective energy of nematic elastomer membranes. A m is the largest principal stretch 
and 5 is the areal stretch. The energy is zero in the region marked C. 


Then, the effective energy of the nematic membrane Wff = W mem is given by 


t < 


,1/3 


hm | i 


,1/2 


5{Ff 


-3 


,i/3 ( *m(F) 2 2 \ 3 

V r A m(F) J 

„!/3 ( ^m(F) 2 + J(F) 




^ 1 

+ 


A m{F) 2 5(F) 2 


if (Xm,5)eC, 
if (Am, 5) e M, 

if (Xm,5) G W, 

— 3 if (Am ,5) € S. 


(5.2) 


Here, 0 < X m (F) < A m(F) are the singular values of F (i.e., the eigenvalues of (F T F) 1 / 2 ), 
5(F) = \ m (F)X M (F) = VdetF T F and 


C := {(X M ,5):X 2 M >5,X M <r 1 / 3 ,5<r 1 / 6 }, (5.3) 

W := {(X M ,5):X M >r 1 / 3 ,5<X]i I 2 }, (5.4) 

M := {(Xm, 5) : 5 > r 1/6 ,r~ 1/2 X 2 M < 5 < X 2 m }, (5.5) 

5 := {(X M ,5):X 1 J I 2 <5<r- 1 / 2 X 2 M }. (5.6) 


Remark 5.2 Some care needs to be taken when dealing with extended real-valued quasiconvex func¬ 
tions. Indeed, the fact that a function f : M 3x2 —> M U {+oo} is quasiconvex (according to the 
defintion of Section of this paper) does not imply that the associated functional f(S7'y')dx' 
is sequentially weak* lower semicontinuous on W 1)0a (uj, M 3 ) m- In the current situation, thanks 
to Remark \f-4l the relaxed energy density has polynomial growth and therefore weak lower sem- 
incontinuity is true for the relaxed functional. Alternatively, we refer the interested readers to Ball 
and James m where a more restrictive definition of quasiconvexity for extended real value func¬ 
tions is presented. This definition guarantees weak lower semicontinuity of functionals associated 
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Figure 2: (a): Idea of the proof of Lemma 5.5. ( b ): level curves of W 2 n in the space The 


web version of this article contains the above plot figures in color. 


to extended real value integrand functions. It is an easy computation to show that both the ap¬ 
proach pursued in what follows and the relaxation technique based on the alternative definition of 
quasiconvex envelope give the same results for the functionals considered in this paper. 


Proof of Theorem 5.1 . Recall that the quasiconvex envelope of an extended value function is not in 


general bounded from above by the rank-one convex envelope. However, we show that this bound is 
true for W 2 n- By Remark 4.4 is a finite-valued, a quasiconvex function and W^ T) = (H 7 |£,) l?c . 

Therefore, if we substitute / = IF^d ' n (2-5) we obtain 


(VF^r < ( W? D ) qc < (W« c D y 


Then, by (2.4) we conclude 

wFd = (w pc D ) pc < (w? D f c < (w&r < (w-TdY 0 < w.f c D 

and recover the classical inequality 

{w 2D ) pc < ( w 2D y c < c w 2D y c . 


(5.7) 


(5.8) 


(5.9) 


We show in Lemma 


Wf c n < W mem 


5.3 


and Lemma |5.4 that W mem < W pr D . We show in Lemma 


5.5 


that 


Combining these with (5.9), 

y^rmem 

and the result follows. 


< W pc D < W.y D < IF l c D < IF" 


(5.10) 


5.2 Step 1: A formula for W 2 d 


Lemma 5.3 For W 2 d defined in (5.1), 


min ipi(\ M (F),5(F)), if rankF = 2 
W 2D (F) = { '6{L-,3} 

+00 otherwise 


(5.11) 
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where 


with 


min cpi(\M,$) 
*e{ l.3} 


<Pi(\m,5) if Am<5 >r 1/2 ,S < X 2 M 
< 9 ? 2 (Am, 6 ) if \ M S <r~ 1/2 ,6 < X 2 M 

^s(A M ,S) ifX M S€ {r- l/ \ r 1/2 ),5 < \\ 


^1 (A m, <5) 
<^2(A Ml <5) 

9^3(AM, A) 


— -f /’1/3 


— -f /’1/3 


3* , , 1 

r + Al / + ^j 


1 


— / t-1/3 


2 ^ r 


+oo 


M+ + r<5 2 J 
<^ 2 o Am 

X 2 + 

A M 


}■ 

*}• 


V 2 <6 


-3 


if \ M d £ (r 1 / 2 ,r 1/2 ), 


otherwise . 


(5.12) 


(5.13) 

(5.14) 

(5.15) 


Proof. The proof is an explicit calculation. To begin, if rank I 7 7 ^ 2, then det(T|c) = 0 for every 
c G M 3 . This implies W e (F\c,n) = +00 for every c G M 3 . Then W- 2 o(F) = +00. Thus for the 
remainder of this section, we restrict our attention to the case that rank I 7 = 2. 

Let fp (c) := Wo(F\c,n) and gp(c) := c T adj F — 1. Then, inf cgR 3 W e (F |c, n) is equivalent to 
the optimization 


j£» {*>) : »*(“) = °} • 


Here fp is a convex, differentiable function and gp is an affine equality constraint. It follows 
that co is a global minimizer of this optimization if and only if there exists a A G M such that 
Vfp n ( c o) + XVgp(co) = 0 (see for instance ( 8 ], Section 5.5.3). Solving this equation, we obtain 


co = 


£ adj F 
\£ 1 / 2 adj F \ 2 


for £ defined in (3.2). 


Let W(F,n ) := Wo(F\co,n). Since cq is a global minimizer for the constrained optimization 


above, W{F,n) = inf cgR 3 W e (F\c,n). Then from (5.1), it follows that inf ng §2 W(F,n) = W 2 d{F). 
For this optimization, we simplify the analysis through a change of variables. We write F = QDR 
for Q G 50(3), R G 0(2) and D a diagonal matrix as in (2.2) with Xm > A r 
values. We can say A m > 0 since rank I 7 = 2. 
constraint via m 2 = 1 — m 2 — mi 2 • Then by direct substitution, 


> 0 as the singular 
Additionally, we set n = Qm, and impose the § 2 


W(QDR, Qm) = g I r 1 / 3 


7A M — £2 (A M - A m) + 


1 


c(7- 1 )AmA^_ 


- 3 > =: ip(\M , A m , 7 , 6 ), 


where 


7 = 6 + 6 , 6 ("h) = 1 - antf, i = l,2, a = 


r — 1 
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Here a E [0,1) since r > 1. Further, we let <5 = A M^m, and set 

p(X M ,S, 7 , 6 ) : = <?(Am,<V a m,7;6) 


= a l r i/3 


7^m — £2 I A^- - ( -— I + 


Am/ / r( 7 -l)J 2 


- 3 


(5.16) 


Note that the constraint A m > A m > 0 implies X 2 M > 8 > 0. 

W is dependent on only four constrained variables. Consider the closed set 

B := {( 7 , 6 ) : 7 - 6 < 1, 6 G [1 - a, 1] and 7 E [2 - a, 2]} . 

£> combined with the constraint A|^ > 5 > 0 give the admissible set for p. Hence, we have 

W 2 D (F) = inf W(F, n) = inf \p(X M (F), 5(F), 7, 6) : A| f > 5 > 0,( 7 ,£ 2 ) G b\ . 

neS 2 7,$2 t j 


Observe that 


inf {p{\ M , <5, 7 , 6 ) : A m>5> 0 , ( 7 ,£ 2 ) E H} = p( X M , S, 7 ,1) =: p 0 {X M , S, 7 ) 

£2 


by (5.16) since the second term in the brackets is non-positive. Then, 

),S{F), 7 ) : X 2 m > 6 > 0,7 E [2 - a, 2] |, 


W 2 D (F) = inf 

7 


where 


<A)(Am, 7 ) = 7 ^ < ?’ 1/3 


(7 — l)A^f + ( -—+ 


A mJ r( 7 — 1)<5 2 


- 3 


is a continuous function on this constrained set (which is moreover bounded in 7 ). It is also 
differentiable for 7 in the open domain (2 — a, 2). It follows that the inhmum is attained. Further, 
7 minimizes po only if it is on the boundary, i.e 7 = 2 — a or 7 = 2 , or it is a critical point, i.e. 
^ 7 ^ 0 ( 7 ) = 0 and 7 E (2 — a, 2). 

We proceed case by case. Letting 7 = 2 — a, observe < po(Xm,8,2 — a) = <^i(A m,S) in (5.13). 
For the other boundary 7 = 2, we obtain pq(Xm,S,2) = p 2 (Xm,S) in (5.14). Finally, in computing 
the critical point d^p 0 ( 7 ) = 0 , we obtain 


7 = 


r 1/2 
A m$ 


+ 1 . 


Direct substitution pq(Xm,S,j) yields the equation given for the finite portion of tps in (5.15). 
Recalling that 7 must lie in the domain (2 — a, 2), we set P 3 = +00 if 7 does not lie in this set. 
This is the full result in (5.11). 

To complete the proof we compute the minimum in Eq. (5.11) thus yielding (5.12). First, 
observe that 


- P 3 = w r 1/3 (r 1 / 2 X M ~ v 


1\2 


Pi - P3 


B 1/3 

<^2 - ¥>3 = 2 r 


A M — 


--1/2v2 


(J / “ 


> 0 = 
> 0 


<P3 < <Pl, 

■ <P3 < <^ 2 , 
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thus proving (5.12) in the region A m5 £ ( r ^ 2 , r 1 / 2 ). To complete the computation in the remaining 
regions observe that 


yielding <pi < ip 2 if S > and ^2 <p\ if 5 < . Therefore we have 

min (pi(\M,S) = ipi if A M 5 > r 1/2 ,5 < \ 2 M 


1\ /I 


1 


and 


min 5) = (f 2 if A M 5 < r 1/2 , 5 < \ 2 M 

1j-j3 


as required. 


5.3 Step 2: Upper bound or W mem < 


Lemma 5.4 Let W r ‘ 


be as in (5.2) and W 2 D as in (5.1). Then, for each F e M 3x2 ; 
W mem {F) < W^(F). 


(5.17) 


Proof. We prove this in two parts. In Part 1, we prove that VP mem is polyconvex and in Part 2 we 
prove that W mem < W 2 D- The result follows. 

Part 1. We now show that l'F mem is polyconvex. First, observe from (5.2) that there exists a 
function if> : Ml —y M (here by M + we denote the set of all non-negative real numbers) such that 

W mem (F) = V>( Xm(F ), 5(F)). (5.18) 


It also follows by verification (also see Proposition 2 of DeSimone and Dolzmann |14| ) that if is 
convex and if) is non-decreasing in each argument (i.e., if(s,t) in nondecreasing in s for fixed t and 
nondecreasing in t for fixed s. We then notice that A m(F) = sup me §i \Fm\ is convex in F. Further, 
5(F) = | adj .F| is convex in adj F. Since the composition of convex function with a non-decreasing 
and convex function results in a convex function, we conclude that there exists a convex function 
g : M 3x2 X M 3 —> M such that 

^(X m (F), 5(F)) =g(F, adj F). 


Combining with (5.18), we conclude that 


W mem (F) = g(F, ad) F) 


(5.19) 


for convex g. By definition of polyconvexity, W mem is polyconvex. 

Part 2. We now show that W mem < IF 2 d■ We show by explicit calculation in the Appendix A that 


It follows that 


from (5.11). 


W mem (F) < ipi(\ m (F), 5(F)), i = 1,... 3. 
r™(F)< min <pi(\ m (F),5(F)) = W 2D (F) 

ie{ 1,...,3} 
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5.4 Step 3: Lower bound or < W mem 

Lemma 5.5 Let W mem be as in (5.2) and W 2 D as in (5.1). Then, for each F e 


3x2 


W^(F) < W mem (F). (5.20) 

The proof makes repeated use of lamination. We collect the calculations in the following proposition. 


Proposition 5.6 Let q, d G M with q> 0, q 2 > d and define 

K. := {GeR 3 * 2 : X M (G) = q,8(G)=d}, 

M w := {GeR 3 * 2 :\ M {G) = q, 5(G) e [0,d]}, 

M s := {GeR 3 * 2 :\ M (G)e[d 1/2 ,q], 8(G) = d}. 

Then, for d > 0 M w c IC^ and for d > 0 M s C IC^l. 

Proof. We begin with M w . Let G e M w with A m(G) = q, 8(G) = <5 £ [0,4]. Using the polar 
decomposition theorem, we can take 


G = 



Define 


G ± = 


q 0 

0 ±f 

0 0 


e=- 1 +- 


Note that G Y± e 1C, 6 e [0,1] since 8 < d, rank (G + — G ) = 1 and G = 9G + + (1 — 9)G . Therefore, 
M w c/C«. 

The proof of M s C K)C> is similar (also see [Q, Theorem 3.1]). Again, using the polar de¬ 
composition theorem, we can take G e M s as a diagonal matrix. First, let us assume c / 0 and 
y/d < c< q which corresponds to c > d/c and define 



(5.21) 


Note 5(G ± ) = d. Further, the eigenvalues of (G ± ) T G ± are 


So the choice 


^le+^+c 2 



^[e+^+c 2 ] -d 2 . 


d 


d 2 


1 




d 


2 , 


q 


c 2 J 


2 2 
q — c 


> 0 


q* c q* 

makes Am(G ± ) = q. Therefore, G^ e 1C. Further, G = \G + + \G~ and rank ( G + — G~) < 1. 
For the case G e M s such that G = 0, replace the diagonal entries in (5.21) with 0 and repeat the 
argument. Therefore, M s C 1C^ 1 \ ■ 
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Proof of Lemma 5.5. We show that 


W^ C D (F) < W mem (F) 


(5.22) 


region by region. In the region S, note W 2 o = W mem and the result follows. 

Now, let FeW with q = A m(F) > ?’ 1//3 and d = 6(F) < g 1 / 2 . This corresponds to the point A 
in Figure [2} Let 

JC = {G G M 3x2 : A m(G) = q, 5(G) = q 1/2 }. 

This set corresponds to the point B in Figure^. By Proposition (jTHJ we have F G . Therefore, 
there exists A G [0, 1] and G\ , G 2 G AC with rank(Gi — G 2 ) < 1 such that 


W^d(F) < 
< 


A WZMGi) + (1 - A )W r 2 c D (G 2 ) 


A W 2D (G 1 ) + (1 - X)W 2D (G 2 ) 


ff 
2 L 


A/3 


r 




Above, the first two inequalities follow from the fact is rank-one convex and Wff < W 2 d- 
The following two equalities is by explicit verification of the formula. 

Now, let F G M with d = 5(F) > r 1 / 6 and q = A m(F) G [d 1 / 2 , r 1 / 4 ^ 1 / 2 ]. This corresponds to 
the point C in Figure^. Let 

K. = {G G M 3x2 : A m (G) = r F 4 d 1/2 , 5(G) = d}. 

This set corresponds to the point D in Figure^. Therefore, by Proposition 
Therefore, arguing as before, Wffj(F) < W mem (F). 

Finally, let F G £ with q = A m(F) < r 1 / 3 and d = 5(F) < min{g 2 ,r 1 / 6 }. This corresponds to 
the point P in Figure [2^,. Let 

K = {G G M 3x2 : A m (G) = r 1 / 3 , 5(G) = d) 


5.6, we have F G 


and 


AC = {G G M 3x2 : A m(G) = r 1/3 , 5(G) = r 1 / 6 }. 


These sets correspond to the points Q a nd Z in Figure [2j-(a) respectively. From Proposition 5.6 


F G Further, again by Proposition |5.6[ K. C AC^\ In othei words, F G IC^ 2 \ We can again 

argue as above to show that Wff->(F) < 0 = W mem (F) as required. 


6 Characterization of fine-scale features 

The energy density W 2 £> is not quasiconvex. Thus a membrane with this energy density is able to 
relax its energy to that of Wf c D through the introduction of fine-scale features. In this section, we 
characterize these features. Briefly, we show that the features in region A4 are essentially planar 
involving oscillations of the director (i.e., no wrinkling) while those in W are necessarily wrinkles 
(i.e., uniform director). Further, we show that there are no fine-scale features in region S. 

To characterize the fine-scale features, we consider the two-dimensional energy 

ho(y) = [ W 2D (Vy)dx' ( 6 . 1 ) 
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subject to affine boundary conditions, i.e the space of deformations Ap := {y E H 1 ^,!! 3 ): y — 
Fx' E i^Q(cu,M 3 )} with F E M 3x2 . It is known (cf. Lemma 3.1(h) and Lemma 6.2, [10]) that there 
exists weakly converging minimizing sequences that satisfy 

yj Fx 1 in H 1 (w,M 3 ) with I 2 D(yj) —>■ inf I 2 d = \w\W mem (F) as j —> 00 . (6.2) 


Let v x be any H 1 gradient Young measure generated by such a sequence. Since {yj} is a minimizing 
sequence for I 2 D, it is also a minimizing sequence for the relaxation f W mem (V'y)dx'. Further, 
since W mem is non-negative and bounded as in (4.9), it follows from Theorem 1.3 of Kinderlehrer 
and Pedregal [20] that 


/(Vyj) (is x , f) in L 1 (A) for any / E { g E C(M 3x2 ): 


sup 

GeR 3x2 


\g(G)\ 

|G| 2 + 1 


< +00 


(6.3) 


and for every measurable A C uj whenever the sequence {/(Vj/j)} converges. As an immediate 
consequence, we obtain the identities 


(Kr, id) = F, 


{v x , W mem ) = W mem (F) a.e. x E u. 


(6.4) 


Now, since W 2 D is a normal integrand, the fundamental theorem of Young measures gives an 
inequality, lim inf hoiyj) > / (z^, W 2 o)dx (cf. Definition 6.27 and Theorem 8.6, [18]). Thus, 


M W mem (F) = lim I 2D ( yj ) > / (; v x ,W 2D )dx > / ( v Xl W mem )dx = \co\W mem (F) 




where we use the fact that W 2 d > W mem . It follows \uj\W mern (F) = f UJ (^ x ,W 2 n)dx. Again using 
the fact that IT 2 D > W mem and (6.4h we conclude 


{v x , W 2D ) = W mem (F) a.e. x E u. 


(6.5) 


By the localizing properties of Ff 1 gradient Young measures (cf. Theorem 2.3 of [19]). we conclude 
that the fine-scale features which arise from minimizing sequences of W 2 d are described by the 
homogenous H 1 gradient Young measures which admit the identities, 


id) = F, (v, W 2D ) = W mem {F). 


( 6 . 6 ) 


We present a characterization of this in the following theorem. 


Theorem 6.1 Let r > 1, F E R 3x2 and let 

Mp := {v E M qc , < 1 /,id) = F, ( 1 ,,W 2D ) = W mem (F )} 


(6.7) 


be the set of homogenous Ft 1 gradient Young measures that satisfy (6.6). Then, there exists u E Mp. 
Further the following is true. 

1. (The region M.) Suppose (A m(F), 5(F)) E M.. Set 6 := 5(F). Let the singular value decom¬ 
position (cf. (2.2)) of F be given by F = QDpR with Q E SO(3),R E 0(2). 

If v E Mp, then 

supp D d fC s := {G E M 3x2 : G = QQD~ S R, Q E 50(3), R E 0(2), det(R)Qf 3 = det (R)f 3 } 

( 6 . 8 ) 
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where E § 2 is orthogonal to the plane of the reference configuration of the membrane and 


4/4 

D s = 5 1/2 ( 0 


,-1/4 


(6.9) 


2. (The region W) Suppose (Xm(F), 5(F)) E W. Set Xm = Xm(F). Further, set cm E § 2 and 
fM £ S 1 to be the unique pair (up to a change in sign) of vectors which satisfy F$m = Xm&m- 

If v E Mp, then 


supp u C K~ x := {G E 


e x2 : (X m , 6)(G) = Gf M = Wm}- 


3. (The region C) Suppose (A m(F),S(F)) E C. 

If v E Mp, then 

supp v C /Co := {G E M 3x2 : (A M ,5)(G) E L with 5/X M = r~ 1/e }. 

4■ (The regionS) Suppose (Xm(F), 5(F)) E S. If v E Mp, then v is a Dirac mass. So suppz^ = 

m- 


Theorem 6.1 has striking physical implications. First consider Part [l] corresponding to region 
Ml and consider the particular case when Q = I. Consider any G E supp v and its characterization 
in (6.8). Since det(i?)Q /3 = det(i?)/ 3 , it follows QD$Rv ■ fs = 0 for each v E M 2 . In other words, 


QDgR maps M 2 to M 2 . Thus, all the oscillations are in the plane. Further, for such matrices G, 

W 2D (G) = W 0 (G\c) = W e (G\c, n) 


for c = (0,0,5 1 ) T and n ■ fo — 0 . 
(Xm(G), 5(G)) E S (see Lemma 


6.2 


The first of these identities follows from the fact that 
below) and Lemma 


5.3 


Now consider Part 


Using an argument as 


while the second follows from the 
fact that the largest principal value of (G\c) is A m(G). Importantly, the director is always in the 
plane. In summary, the director oscillates in the plane and oscillations create no out of plane 
deformation. The case Q ^ I is similar except the plane is oriented by the rotation Q. Thus, the 
fine-scale features in /VQs limited to in-plane oscillations of the director. 

/ f 

First consider the case when e,v/ = ( 

before, for any G E supp u, 

W 2D (G) = Wq(G\c) = W e (G\c , n) 

for c • eM = 0, |c| = X M and n = cm- In other words, the director n is fixed with an in-plane 
direction eM- Further, notice G is necessarily of the form 

0 

{ -l/2 

w 

0 



in the eM ~ 


f-L \ 

q 1 J — fs frame for Q E 50(3) that satisfies QeM = eM- In other words, the 

membrane is uniformly deformed and the fine features are related to rotations about a fixed axis 
eM- In other words, oscillations represent wrinkling and these oscillations are always perpendicular 
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Im 

0 


is similar except a uniform rotation orients 


Im 

0 


to 


to eM- The general case eM 
eM- 

Part [3] says that region £ involves only the spontaneously deformed states while Part [4] says 
that there are no fine-scale features in S. 

We now turn to the proofs of the theorems. They rely on the following lemmas. 

Part [7J Then any v E Mp 


Lemma 6.2 Let F E 

satisfies 


3x2 


and 6 satisfy the hypotheses in Theorem 


6.1 


Lemma 6.3 Let F E 

P E Mp satisfies 


supp P c{Ge M 3x2 : (A m,S)(G) = {r l / A 8 1/2 ,5)}. (6.10) 

M 3x2 and Am satisfy the hypotheses of Theorem 6A_ Part Then any 


supp v C {G E 


1)3x2 


:(Am,5)(G) = (Am,A^ 2 )} 


( 6 . 11 ) 


Proof of Lemma [6J% Recall from Section [ 5 ] that we may write W mem = ip o (A m,8) and Wfin = 
p o (A m,$) where ip (p) : 7Z —>■ R (R U {+00}) respectively for 7Z = {( s,t ) E M 2 : s 2 > t,t > 
0}. Recall also that ip is a convex, and it is non-decreasing in each argument. Also, ip < p. 
Finally, (Am,< 5) : M 3x2 —y 7Z are quasiconvex functions bounded quadratically. Therefore, for every 
homogenous H 1 gradient Young measure with (u, id) = F. 

W mem (F) = iPo(\ M ,6)((v,id)) 

< Am), (v,6)) 

< {v,ip° (Am, 5)) 

<(iy,po(X m ,5)) = (v,W 2 d)- (6-12) 

Here, the first inequality follows from the Jensen’s inequality satisfied by homogenous H 1 gradient 
Young measures since (Am, <5) are quasiconvex with the appropriate growth and ip is non-decreasing 
in each argument. The second inequality follows from the convexity of ip , and the third follows 
since ip < p. 

Now, for any v E Mp, each inequality in (6.12) is an equality. This restricts the support of u. 
To deduce this restriction, suppose that the point (A m(F),6(F)) E M corresponds to point C in 


Figure 2a 


Consider the first inequality. By quasiconvexity and growth conditions, (u,X m) > A m(F) and 
{p, 5} > 6(F). In the Am _ 5 space in Figure 2a these inequalities imply the point ((u, Am), (P, 8}) 


cannot be to the left or below point C. Further, every point to the right and above the point C has 
higher ip (cf. Figure [I]) except the line between and including the points C and D. Hence, 


((y Am), (y ^)) £ CD. 


(6.13) 


Next, consider the last inequality. Since p = ip only on5ll{(s,f) E C : t/s = r 1 / 6 } =: S' (see 
Figure^), we conclude 


supp v C \ G E 


{g E M 3x2 : (X m ,5)(G) £5'}. 


(6.14) 


It remains to consider the middle inequality in (6.12). We do this in Proposition 6.4 below. If 
the middle inequality is an equality, we show in the proposition the support of P satisfies (6.10). 
This completes the proof. 
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Proposition 6.4 Let F and 6 be as in the Theorem 6.1 Parf[7| If u satisfies ( 6.13), (6. 14) and 

ip((0, Am), (P, 5)) = (P,ip o (Am, 5)), (6.15) 


then the support of v satisfies (6.10) in Lemma 6.2. 


Proof. Set A + = {G: (A M , S)(G) € S n {5 > 5}} and 


0 + = 


IA+ 


dP(G). 


If 9 + = 1, then by the polyconvexity of 5, (v,5) > 5 contradicting (6.13). Now consider the case 
1 > 9 + > 0. Set 


A m e+ 


4+ 


A M (G)dP(G), 


6 + := 


9+ 


'A+ 


5(G)do(G), 


A" - 1 

a m - i _ e+ 


5~ : = 


2 \A+ 


A M (G)dP(G), 


1 - 9 + 


2 \A+ 


5(G)d0(G). 


Clearly, 5 + > 5 and 


0 +A+ +(1-0+)Am = (Mm>, 

9 + 5 + + (1 - 9 + )5~ = (p,5). 


From the equality in (6.16) ,6~ <6. Further, notice from the convexity of ip that 


^(A^,6 + ) < ± I if (^m(G), 5(G) J dP(G), 




V»(A M ,S ) ^ 


1 


1-9+ 


2 \ 4 + 


ip 


dP(G) 


(6.16) 

(6.17) 

(6.18) 


Now, in the A m~5 space shown in Figure [2a| the definitions above imply that the point (A^, d + ) 
is a point above the line CD while (A Jj,5~) is below the line CD such that the line joining these 
points intersect CD. It is easy to verify by explicitly computing the derivative along such lines (or 
by inspecting Figure [l]), that if is strictly convex in such segments. Therefore, 


^((^ Am), (0,5)) = (^ + Am + (! - 9 + )\- m ,9 + 5+ + (1 - 9+)6~) 

< 9 + fj( A+, 6+) + (1 - 9 + )if{S~) 

<(P,ip° (A M,S)). 


(6.19) 


The last inequality follows from (6.17). However, this contradicts the assumption (6.15). 
Therefore, 9 + = 0, and 


supp P c {(Am,<5)(G) e S' : 5(G) < j}, 


( 6 . 20 ) 


which is the comp limen t of A + in the set given i n (|6. 14 ). 

Finally, given (6.20) and since 5 = (0,5) ( see 6.13), it follows that supp 0 C {( Am, 5 )(G) E S' : 
5(G) = 5}. But this is just a single point in the Am — 5 space, and it’s given by (6.10). Thus, we 
conclude the proposition. 
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The proof of Lemma 6.3 is very similar and omitted. 


Proof of Theorem 6.1. Existence of a v E Mp follows from the construction in Section 5.4 
Part[7| For any F with (A m(F), 5(F)) E M. and for any v E Mp, the support of u satisfies (6.10) 


by Lemma 6.2 Note that (Xm,5)(G) = (r 1 / 4 ^ 1 / 2 , 5) is equivalent to stating that the principal 


values of G are r 1 / 4 ^ 1 / 2 and r V 4 ^ 1 / 2 . Therefore, by the singular value decomposition theorem 


(2.2), it follows that 

supp v C {G E 


p3x2 


: G = QD S R, Q E 50(3), R E 0(2)} =: K. 


supp 


( 6 . 21 ) 


for Dg is given in (6.9). 

Now, for any D E M 3x2 , Q E 50(3), R E 0(2), it is an easy calculation to find that adj (QDR) = 
det(P)<5adj D. Further for D of the form (|2.3[), adj D = A 1 A 2 / 3 . Further, the adjugate is a minor 


and therefore ( 1 /, adj) = adj (( 1 /, id)) = adjP. Recalling the support (6.21) of v, we conclude 


det (R)Qh = i(^,adj) = i 

0 0 


adj Gdu(G) = \ 
3x2 Jk-s 


det (R(G))Q(G)f 3 dn(G). 


Note that det(R)Qfs E § 2 , and det (R(G))Q(G)f 3 E § 2 for each G E K supp . In other words, the 
equation above states that an average of a distribution on § 2 yields an element of § 2 . However, 
since each element of § 2 is an extreme point, it means that the distribution is concentrated at a 
single point on § 2 . That is, if we let Qo(G) = Q T Q(G), then det(P )/3 = det(R(G))Qo(G)f 3 . The 
result follows. 

Part [^1 For any F with (Xm(F), 5(F)) E W and for any v E Mp, it follows from the definition of 

e-M, fM that 


GfMdv(G) = A m&m- 


( 6 . 22 ) 


So, 


\Gf M \dv(G) > 


Gf M du(G) 


= | P/m | = Am- 


However, from Lemma 6.3, we see that max eg §i \Ge\ = Xm(G) = A m for each G E supp v. There¬ 


fore, | Gfu | < A m f° r each G E supp v. We conclude th at \GJ m\ = A m for each G E supp v. 
Setting G'/m = A m&(G) for e(G) E § 2 and substituting in (6.22), we conclude that e(G) = for 
for each G E supp u. The result follows. 

Parf[s| For any F with (A m(F), 5(F)) £ A the result follows from the fact that W- 2 D is non-negative 
and W 2 d(G) = 0 if and only if G E /Co- 

Part[7} Finally, let F E M 3x2 such that (Xm(F), 5(F)) E S, 9 E Mp. Recall W 2 D = T 0 (A m,5) 
and ip is strictly convex in 5. Thus, 


supp v C {G E 


p3x2 . 


: ( X M ,5)(G) = (A m ,S)(F)}. 


This is actually equivalent to the set (6.23) given in Proposition 6.5 below since X m = 5 /Xm■ The 
result follows from the proposition. ■ 


Proposition 6.5 Let F E M 3x2 such that the singular values satisfy the strict inequality 

Xm(F) > X m (F) > 0. 
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Suppose is is a probability measure on the space o/M 3x2 matrices such that (is, id) = F and 

suppi/ C {G G M 3x2 : (Am, A m )(G) = (A m , A m )(F)}. (6.23) 

Then is (up to a set of measure zero) is a Dirac mass at F. 

Proof. To begin, set (Xm,Xu) = (Xm(F), X m (F)). We let {ei,e2} C M 3 and {/i,/2} C M 2 be sets 
of orthonormal vectors such that 


F = X M ei <8> /i + A m e 2 <g> fi- 


(6.24) 


Let <pj ] (G) := \Gfi\ 2 . This is a convex function. Therefore, by Jensen’s inequality and given 
(is, id) = F with F satisfying (6.24), 


> TpiF) = Am- 


(6.25) 


Conversley, applying a similar change of variables (6.24) to the G G suppz', we see 

= J (KAmci <8> fi + X m e 2 <8> / 2 )/i | 2 ) (G)dis(G) 

' Am cos (9(G)) 2 + X 2 m sin(0(G')) 2 ) dis(G) 


f = X 2 m if is({G G M 3x2 : sin(0(G)) 0}) = 0 

I < X 2 m otherwise, 

since by assumption Am > A m. He re, cos 6 denotes the direction cosine between fi and f\. Com¬ 
bining this observation with (6.25), we deduce (up to a set of measure zero), sin(0(G)) = 0. This 
implies (up to a change in sign) f \ = f\ in measure. Since f\ and f 2 are orthogonal, it follows that 
(up to a change in sign) f 2 = f 2 in measure. 

We repeat this argument substituting with the convex function ip Sl (G) = \G T e \\ 2 . It follows 
that (up to a change in sign) e\ = e\ and e 2 = e 2 in measure. The fact that (is, id) = F ensures 
the eigenvectors are Hxed and not oscillating in sign with some non-zero measure. The conclusion 
follows. ■ 


7 State of stress and connection to tension field theory 


In this section, we seek to understand the state of stress in the membrane. 

Formally, consider an incompressible energy density W%d of the form in (3.4) and assume Wq 
is C 1 differentiable. The Piola-Kirchhoff and the Cauchy stress are defined as 


P(F) = X7 f W 0 (F) - p( adj F) t a(F) = (X7 f W 0 (F))F t - pi 


(7.1) 


where p is the indeterminate pressure (Lagrange multiplier to enforce incompressibility) and I is 
identity. We find p by requiring the tractions to be zero on faces of the membrane. Alternately, 
recall that we obtain the membrane energy density W 2F by writing F = (F\c) and minimizing with 
c (when F is full rank). The minimizer cp satisfies 

S7 c W 0 (F\cp) - p(adjF) = 0, cp ■ adj F = det F = 1 p = X7 c W 0 (F\cp) ■ cp. (7.2) 
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Above, V c denotes derivative with respect to the third column of the deformation gradient. Sub¬ 
stituting this back in (7.1) and writing VfM'o = (VpWo|V c Wo) we obtain a characterization of 
the state of stress in the membrane. 


PoMF) := P(F\cp) = (Vj,Wo|0) = (VpW 2D \0), 

<t 2 d(F) := a{F\cp) = (V pW 2D )F T . (7.3) 

Notice that these depend only on Win- 

However, the effective energy of the membrane is not W- 2 D but its relaxation. In other words, 
energy minimization with the integral of W 2 D can lead to fine-scale oscillations, and thus the 
stress may also oscillate on a fine scale. Therefore, we need to understand the overall of effective 
stress. Ball et al. [2] have shown that if / : M nxm —> M is differentiable and satisfies certain 
growth conditions, then f qc is a C 1 function. Moreover, V/ 9C can be written in terms of V/ and a 
homogeneous W l ' p gradient Young measure u generated by minimizing sequences of j ^ f(Vy)dx , 
i.e. 

Vf qc = J Vfdu. (7.4) 

Unfortunately W 2 D is an extended function (equal to +00 when rank F < 2), and the analogous 
result is unknown. However, our resulting effective energy = \Y mem i s finite everywhere and 
is differentiable except on a boundary. So have the following characterization of the stress. 


Theorem 7.1 Let r > 1, let V C 
let F E 


5x2 


0 < t < s 2 }, and 


l 2 be the open set V := {(s,t) E 

such that (Xm(F), 5(F)) E V. If up is a homogenous H 1 gradient Young measure 
generated by minimizing sequences for the energy I 2 D in the space Ap (see (6.1)) with support in 
V. then 


V pW mem = / V d W 2D dup(G), 


(V pW mem )F T = I (X7 qW 2 d)G 1 dup(G). 


(7.5) 

(7.6) 


Further, the Cauchy stress a mem (F) := (V pW mem )F T has the following explicit characterization. 


(0 


a 


mem _ 1/3 < 


- i ) Id2 


- 1/2 <52 




1 


r A M 


ei <g> ei 


*/( AM,5)e£nD, 
if(\ M ,S) eMnv, 

if(X M ,S) e wnp, 


(7.7) 


A 2 


— - G ) e\ ® e\ + ( -2-yy ) e 2 <8> e 2 if (Am, 6) E S n V. 


1 


<5 2 


1 


A m 3 2 


Before we prove the theorem, we make a few comments on the physical implications. First, the 
membrane is always in a state of plane stress in the tangent plane. Second, the principal stresses 
(the eigenvalues) are always non-negative. Therefore, the membrane can not sustain compressive 
stress. Further, the stress is zero in region £, uniaxial tension in W, equi-biaxial tension in JY[ and 
biaxial tension in S. The different regimes are shown in Figure [3j 
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Figure 3: The effective stress of nematic elastomer membranes. 


To understand this further, consider the special case r = 1 when this theory reduces to that of 
neo-Hookean elastic membrane. The region A4 now disappears and we are left with regions C, W 
and S with zero, uniaxial tension and biaxial tension respectively as in the traditional tension field 
theory [25j [26l [28] . 

Nematic elastomers membranes with r > 1 are characterized by an additional region Ai 
where the state of stress is equi-biaxial tension. This is true even though, the principal stretches 
(A Mi 6 /Xm) can be unequal. In other words, one can have shear strain but no shear stress. This is 
a potentially useful attribute of liquid crystal elastomers in membrane applications. 

We turn now to the proof of Theorem |7.1| 


Proof of Theorem 7.1, Recall from Section [5] that we may write W 2 D = T 0 (Xm,6) and W r 
if o (Am, 5). Now, any G £ M 3x2 has the representation 


G = \ M ei <g> /1 + A m e 2 <8> / 2 , 


(7.8) 


where {ei,e 2 } C M 3 and {/i,/ 2 } C M 2 are orthonormal and A m > X m > 0 are the singular values 
of G. These singular values (and therefore 6 = X m X m) are continuously differentiable with respect 
to G as long as they are distinct, i.e. A m > with 


V^Am = ei <8)/ 1 , V^A m = e 2 <g> / 2 , 

(cf. Corollary 3.5 and Theorem 5.1, |27]). We can use this fact and the representation for (p, if in 
Theorem |5.1| and Lemma |5.3| to conclude that W^n and W mem are continuously differentiable on 
{G : (X m {G),6(G)) £ V}. 

The rest of the proof is by computation and verification. 


Case 1: (Xm{F),5(F)) £ Ai C\V. Set 5 = 6(F). According to Theorem 6.1 Part[lJ suppM^ C /Cj. 


We can now apply the representation (7.8) to F and G £ supp^ to write the identity F = (up, id) 


as 


Xm^i <8> fi + (^J^j e 2 <8> /2 = 6 1/2 J ^r 1/4 ei <8> fi + r 1/4 e 2 ® / 2 j (G)dup(G). (7.9) 
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Another implication of Theorem 6.1 is any G £ supp up can be written as QQD$R where Q £ SO( 3) 


arises from the identity F = QDpR, for some Q £ SO(3 ) and R £ 0(2) such that det(i?)Q / 3 = 
det(R)f 3 . Here, / 3 £ S 2 is orthogonal to the reference configuration of the membrane. Without loss 
of generality, we assume / 3 = (0,0, 1) T . Now for each a = 1, 2 there is a corresponding c a > 0 such 
that e a • 0 Qh ) = c a (QQD- s Rf a ) • (Q/ 3 ) = c a (D~ s Rf a ) • (Q T / 3 ) = c a (det i?/ det R)(D~ s Rf a ) ■ / 3 = 0. 
In other words, the vectors ei and e 2 span the plane perpendicular to Q / 3 for each G £ supp up . 
Moreover Q /3 = ei x e 2 , and therefore e\ and &2 also span this plane. Now, let Rq be a 90 degree 
rotation about / 3 and Qo be a 90 degree rotation about Q / 3 so that Rq fi = / 2,^/2 = — fi and 
Qoei = < 32 , QoC 2 = — e\. Since e\ — &2 span the same plane as e\ — <= 2 , and /1 — f 2 the same plane 
as /1 — / 2 , we have the following relation. If det(.R) = 1, then 


If det(R) = — 1, then 


-R0/1 =/2, Qoei = e 2 if det(R) = 1 , 

Fq fi = fit Qoei = -62 if det(i?) = — 1 . 

R^fi = -f 2 , Qoei = — e 2 , if det(i?) = 1 , 
Roll = /2, Qoei = e 2 , if det (R) = - 1 . 


(7.10) 


(7.11) 


Thus, pre-multiplying and post-multiplying the identity in (7.9) by 6 1 ^ 2 Qo and Ro respectively 
yields the identity 


( fi /2 )^®/ 2 + 


'P/ 2 ' 

Am 


ei (8) fi = 


lie, 


He, 


J/4 Q 0 ei ® RT fl + r- 1 /4Q oe2 ® / 2 ) (G)d^(G) 


r 1(/4 e 2 <S> f 2 + r ^ 4 ei <8> / 1 ) (G)dup(G) 


(7.12) 


by fl7lQ| ) and ( |7T1 ). 

It is easy to verif y ^^ -(r 1 / 4 ^ 1 / 2 , 5) = 0. Thus, combining explicit differentiation evaluated i 
/Cj with the identity (7.12), we observe 

[V d W 2D dup(G) 


m 


= r 


rh n 

VqAm + ) dup(G ) 


t— ei <8> /1 + A M e2 <8 /2 
Am 



(G)dup(G) 

) (G)d^(G) 

fi + r 1/4 e 2 (8) / 2 ) (G)dup(G) 


= /rr 1 / 3 


e 2 <8 /2 


- e 2 <8) / 2 


(7.13) 


Finally, it can be verified explicitly that V pW mem coincides with (7.13). This gives the identity 
(7.5[) for region M (I'D. 
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Similarly, 


>K.j 


V d W 2D G T dup{G) = I f 6 5 ) G T dup{G) 


= /rr 1 / 3 


= /rr 1 / 3 


= /rr 1//3 


IK* 


1 


^2 - p j ( 6l ® ei + e 2 <8> e 2 ) (G)dup(G) 

Jlp ~ J ( e i ® ei + e 2 (g) e 2 )(G)dup(G) 
r ^ 2 ~ p) W2 ’ 


(7.14) 


The fourth equality uses the fact that the basis {ei (G), e 2 (G)} a lways spans the same plane. Finally, 
it can be verified explicitly that V pW mem F T coincides with (7.14) in this region. This gives the 
identities ( |7.6[ ) and ( |7.7[ ) for region M Pi V. 


6.1 


Part[2j supp up 


C 


Am 


Case 2: (A m{F), 8(F)) E WnV. Set A m = A m(F). Following Theorem 

and so any G E supp z^p satisfies (Xm(G),8(G)) = (Am,A^). In addition, for the vectors f\ E S 1 
and e\ E S 2 such that Ffi = e±, G E supp up also satisfies Gf\ = e\. Writing G E supp up as in 
(7.8), we observe using the properties of the set AN , 


Gfi = (Xm^i <S> h + (5/AM)e 2 <8> /2)/i 
= (Amci ® fi + A m e 2 <8> f 2 )fi 

— —_-j /o — 

= Am cos(0)ei + X M x sin(0)e 2 = A M ei. 


(7.15) 


He re, co s(0) denotes the direction cosine from f± to f±. Applying the squared norm to the identities 
(7.15) yields |G/i| 2 = (Am) 2 cos(0) 2 + X~^ sin(0) 2 = X 2 M . Since X 2 M > X ^ in W, we deduce from 


m 


this equation that cos($) = ±1. That is, f\ is up to a change in sign equal to f\. Substituting for 
fi back into (7.15), we find ei = ±ei when fi = ±/i, or alternatively 


ei <8> fi = ei (g) /i V G E supp up. 


(7.16) 


Now, it is easy to verify explicitly ^j(Xm, X\j 2 ) = 0. Thus, combining explicit differentiation 
evaluated in AN with (7.16), 


V d W 2D dup(G) = 


>K~, 


() dup{G) 


= pr 1 ' 3 


= pr 1 / 3 


= 


\dX M ' G 

Am 




5 2 


Am 

r 

Am 

r 


r 

1 








ei 


hj (G)dup(G) 
ei <E> fi(G)dup(G) 

r 

ei (g> fi. 


(7.17) 


Finally, it can be verified explicitly that V^IF mem coincides with (7.17). Therefore, the identity 
(17.51) is satisfied for W. 
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Similarly, 


V d W 2D G T du p (G) = 


>K., 


(£ V < 3 A " ) 67 ’ <i ^ (G) 


= fir 1 / 3 


= fir V 3 


lie. 


—Y - ^ 2 ~ ) ei (8) ei(G)dvp(G) 
r a m< 


— ~ T~ ) I ei ® ei{G)dvp(G) 


Am 


//Ct 


= /ir 


- ()T .i/3 / AK_^ 


Am 


ei <x> ei. 


(7.18) 


For the last equality, recall e \ = ± ei for G G supp vp. Finally, it is easy to ve rify explici tly t hat 
V pW mem F T coincides with (7.18) in this region. Thus, we have the identities (7.6) and (7.7) for 
region W fl D. 

Case 3: (\m(F), 5(F)) G C n V. According to Theorem 6.1 Part [3j suppi^ C /Co- We see that 
V1T 2D = 0 on K o and similarly \7W mem = 0 on C Fl V. The identities (7.5), (7.6) and (7.7) for 
region C n V. 

Case 4-' (^m(F), 5(F)) G S CD. According to Theorem 6.1 Part [4j supp^ is a Dirac mass. 
According to Theorem |5.1[ W 2 d and W mern coincide on SCD. The identities (7.5), (7.6) and (7.7) 
for region S CD. ■ 
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A Appendix A: Proof of Part 2 of Lemma 5.4 


A.i w mem < (p 2 

First, in the region £ of liquid behavior there is nothing to prove. Therefore, referring to Fig. [4j 
we are left with showing that W mem > (p 2 in the light gray region of equations for r 1 / 3 < Am < 
r~ 1 / 2 S~ 1 . Recalling that in this region W mem = %[r 1 / 3 ((\M(F)) 2 r~ 1 + 2(Am(-U))~ 1 — 3], it is 
enough to prove that 

¥+ jm i y (( a «+¥ + ¥) for: a s r ~ i/2A u am > ri/s }' (ai > 

The critical point of X 2 M + 5 2 \J 2 + r~ 1 5~ 2 is attained at 5 2 = Am?’ - 1 ^ 2 - This corresponds to a 
minimum, yielding the following inequality 


Am 2 o 
~ + t— < A^ 
r Am 


+ 


A Mr 1 / 2 


(A.2) 


34 

















r 1 / 3 Am 


Figure 4: Regions where the comparison for VF mem and the functions ip\ (dark gray), if 2 (light 
gray) and ip% (silver/intermediate gray) occurs. The web version of this article contains the above 
plot figures in color. 


which is indeed true for Am > r 1 / 3 . Then, evaluation of X 2 M + 5 2 X h 2 + r 2 along the curve 
5 = r~ 1//2 X~^ does not improve the inequality above. 


A.2 W mem < (p x 

We focus on the interval 5 > r 1//2 A^ corresponding (if again we ignore the region C) to the dark 
gray area in Fig. [4j This set has a non-empty intersection with both the simple-laminate regions 
M, W and the regime of solid behavior S. First of all, notice that if (Am, 5) € S then W mem = W 2 D 
and there is nothing to prove. 

Let us assume r _1//2 \ 2 M < 5 < X 2 M , 5 > r 1//2 A^/. This corresponds to a subset of A4 for which 
we have W mern = ^[r 1 / 3 (25(F)r~ 1 / 2 + ( S(F ))~ 2 ) — 3]. We are left with the inequality 


+ 4<^ + ^ + ^ for r- 1/2 \ 2 M < 6 < \l„S> 


25 


*1/2 5 2 


A M 


which is trivially true. 

Then, let us assume r 1//2 XjJ < 5 < A^ 2 . This is a subset of the region W for which we have 
W mem = f {rV 3 [((A M {F)) 2 r~ l + 2(Am(F)) _ 1 ] - 3}. The inequality 


A M 


2 A 2 , 

_l-< 

r A m r 


+ + for rl/2 A m < s < a m 2 > 


a m 


follows trivially. 


A.3 W mem < 

We now focus on the interval 5~ l r~ 1 / 2 < Am < 5~ l r 1 / 2 corresponding to the silver/intermediate 
gray area in Fig. [4j Notice that if we remove the region C (for which there is nothing to prove), 
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we are left with two disjoint subsets. 

We begin with considering A m > r 1//3 . Since in this region W mem = ^[r 1 ^ ((X M (F)) 2 r~ l + 
2{\m{F))~ 1 — 3], it is enough to show that 


Am 

which is equivalent to 
2 


\ 2 

a m 


+ T- < 

r Am 


inf| (\m + 2 * ) for -^-r 1/2 < X m < -^-r 1/2 , A M > r 1/3 |. (A.3) 

1 ' r V' a m a m 4 


The critical point of A 2 n + 2r~ 1 / 2 A“ 1 is attained at (A m , Am) = (r -1 / 6 ^ 1 / 3 ) E L. Then, we have 
to evaluate A^ + 2r _1 / 2 A“ 1 on the curves of equations A m = A^/r" 1 / 2 and A m = A^/r 1 / 2 . Notice 
that, for A m = r 1 // 2 AA 2 we have that <£>3 = ipi while for A m = r _ 1 // 2 A7 , 2 we have 923 = 992 and from 


discussion of these cases in Paragraph A.2 and A.l, respectively, it therefore follows that (A.3) is 
true. 

To conclude, we have to prove that the inequality VP mem < <^ 3 holds in the remaining subregion 
defined by X 2 M > 5, 5 > r 1 / 6 and 5 < r 1 / 2 A^ f 1 . This subset is contained in the region A4 in which 
case we have W mem = ^[r 1//3 (25(F)r -1 / 2 + (5(F))~ 2 ) — 3]. Therefore, we are left with proving the 
inequality 


25 

^72 


+ h- inf {(4 + Si) for 61,2 - Xm - s> rl/8 " 5 < SS- 


equivalent to 

25 


+ ^2 - inf { ( A m + r i/ 2 A ) for A m < 5 1/2 ,5 > r 1/6 , X m > r 1/2 5 2 }. 


(A.4) 


r 1 / 2 5 2 “ IV m r 1 / 2 X r 

In order to prove the inequality above it is enough to evaluate the function A 2 , + 2r _1 / 2 A“ 1 on 


the boundary of the region defined on the right hand side of (A.4). This yields the following two 
relations 


25 1 / , 

+ < 5 + 


-1/2 52 


7-1/2(51/2 


for 5 E (r 1//6 , r 1//3 ) 


(A.5) 


^ + 4<r-‘i- + | fbr* 6 (r‘/«, r V>), 


(A.6) 


obtai ned b y evaluating A^ + 2r~ 1 / 2 A“ 1 for X m = 5 1 / 2 and A m = r _1 / 2 5 2 respectively. To show 
that (A.5) holds it is convenient to operate the change of variable (r 1 / 4 ,?- 1 / 2 ) B y := 5 3 / 2 and thus 
writing (A.5) as follows 


y 


2 (r 1 / 2 — 2) + 2y — r 1//2 >0 for y E (r 1//4 , r 1//2 ) 


(A.7) 


which can be easily shown to be true Vr > 1. Then, it is immediate to prove (A.6). 
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